APS/123-QED 



O 
O 

(N 

1—5 






o: 
-(— > 

^ : 

O ■ 

• ^H . 

£?■ 

(n: 

>■ 

m ■ 
en ; 

(N . 

^^ 

O' 

o. 

o ; 

•I— I 

Oh. 



X 



Stark effect of the cesium ground state: 
electric tensor polarizability and shift of the clock transition frequency 
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(Dated: February 2, 2008) 

We present a theoretical analysis of the Stark effect in the hyperfine structure of the cesium 
ground-state. We have used third order perturbation theory, including diagonal and off-diagonal 
hyperfine interactions, and have identified terms which were not considered in earlier treatments. 
A numerical evaluation using perturbing levels up to n = 18 yields new values for the tensor 
polarizability a2(6S'i/2) and for the Stark shift of the clock transition frequency in cesium. The 
polarizabilities are in good agreement with experimental values, thereby removing a 40-year-old 
discrepancy. The clock shift value is in excellent agreement with a recent measurement, but in 
contradiction with the commonly accepted value used to correct the black-body shift of primary 
frequency standards. 
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Since its discovery, the Stark effect has been an im- 
portant spectroscopic tool for elucidating atomic struc- 
ture. The Zeeman degeneracy of the nSi/2 ground state 
of alkali atoms cannot be lifted by a static electric field 
because of time reversal invariance. However, the joint 
effect of the hyperfine interaction and the Stark interac- 
tion leads to both F and |M| dependent energy shifts 
which in cesium are 5 and 7 orders of magnitude smaller 
than the shift due to the second order scalar polarizabil- 
ity. While the scalar Stark shift is understood at a level 
of (1-2) • 10~^ HI SIS) there has so far been no satisfac- 
tory theoretical description of the F and \M\ dependent 
alterations of the scalar effect at the levels of 10^^ and 
10~^. In this Letter, we extend previous theoretical mod- 
els which allows us to bridge a long standing gap between 
theory and experiment. 

The effect of a static electric field on the hyperfine 
structure was treated in a comprehensive paper by Angel 
and Sandars 0| , who showed that the second order Stark 

effect of a state I7) = \nLj, F,AI) can be parametrized 

(2) 
in terms of a scalar polarizability, Oq (7) , and a ten- 

sor polarizability, aj (7), where the latter has non-zero 
values for states with L > 1 only. As a consequence, 
the spherically symmetric nSi/2 alkali-atom ground state 
has only a scalar polarizability, so that its magnetic sub- 
levels \F,M) all experience the same Stark shift, inde- 
pendent of F and M . However, an experiment by Haun 
and Zacharias in 1957 lq| showed that a static electric 
field induces a quadratic (in field strength E) shift of 
the |F = 4, M = 0) <-^ |F = 3, M = 0) hyperfine (clock) 
transition frequency (i^-dependent effect). In 1964, Lip- 
worth and Sandars (g demonstrated that a static elec- 
tric field also lifts the Zeeman degeneracy within the 
F = 4 sublevel manifold of the cesium ground state (in- 
dependent effect). An improved measurement was per- 



formed later by Carrico et al. [3 and its extension to five 
stable alkali isotopes was done by Gould et al. [3]. 

Sandars :'9j has shown that the F- and M-dependence 
of the Stark effect can be explained when the pertur- 
bation theory is extended to third order after includ- 
ing the hyperfine interaction. The theoretical expres- 
sions for the tensor polarizabilties 0:2 of ^ were evalu- 
ated numerically in [g and [g under simplifying assump- 
tions. Comparison with the experimental polarizabilities 
showed that the absolute theoretical values were system- 
atically larger for all five alkalis studied in |^ . Our recent 
measurements [T^, llil [l2l of ao of ^'^^Cs confirmed the 
earlier experimental values 0j lal- To help clarify this 
long-standing discrepancy we have reanalyzed the third- 
order calculation of the Stark effect in the alkali hyperfine 
structure. We have identified contributions which were 
not included in earlier calculations and as a result we ob- 
tain theoretical values which are in good agreement with 
experimental results [3, H, llfl HB UJI ■ When applied to 
the Stark shift of the clock transition frequency our cal- 
culations yield a value in good agreement with a recent 
measurement |13J . 

The third order perturbation of the energy Ea of the 
sublevel la) = |6S'i/2,-F, M) is given by 
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where Ep and E^ are the energies of unperturbed states, 
and the perturbation W — H^f -I- Hgt describes the hy- 
perfine and Stark interactions. Of all the terms obtained 
by substituting W into Eq. ^ only those containing the 
product of two matrix elements of Hst and one matrix 
element of Hht give nonzero contributions because of the 
selection rules AL = ±1 for Hgt and AL = 0, ±2 for H^. 
We address first the second term of Eq. |^. The 
diagonal matrix element in front of the sum represents 
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FIG. 1: Contribution of diagonal hyperfine matrix elements 
to the third order Stark effect calculation. The dotted lines 
represent Stark interaction matrix elements while the solid 
lines represent hyperfine matrix elements. 



only the Fermi contact interaction in the ground state, 
while the sum is carried out over Stark interactions only. 
Diagram A of Fig. ^ represents this term in graphical 
form. The hyperfine matrix element (a| H^f |a) is scalar, 
but F dependent, and the sum is similar to the ex- 
pression for the second order polarizability, except for 
the squared energy denominator. As in the theory of 
the second order polarizability |j| the sum can there- 



fore be expressed as (a| Hj;ff 



(fc=0) 



H, 



(fe=2) 



a). The scalar 



contribution, (ajH^ff |a), depends only on the strength 
of the applied electric field, E^, while the second rank 

("2) 

tensor contribution H^,/ depends on its orientation as 
3E^ — E^. The selection rules for tensor operators re- 

quire that (S'1/2 || H^,/ || S'1/2) =0. As a consequence the 
second term of Eq. |^ gives a scalar contribution to the 
energy which depends on _F, but not on Af , and which 
can be parametrized by a third order scalar polarizability 
a[,'^(65i/2,F)as 



r(3)/ 



1 



.,(3)/ 



t.E'^>{F) = ---c.'^'{F)-^\ 



(2) 



This term produces the major contribution to the Stark 
shift of the hyperfine transition frequency ^vqq — v(F — 
4,M = 0)-i/(F = 3,Af = 0). 

We next address the first term of Eq. (QJ and consider 
only diagonal matrix elements of Hj,,. As before, the 
Stark interaction part of the first term has only a rank 
fc = (scalar) contribution. The dipole-dipole and elec- 
tric quadrupole parts of the hyperfine interaction have 
the rotational symmetry of fc = 0, 2 and k = 2 tensors, 
respectively. Together with the scalar Stark interaction, 
the first term thus has scalar and tensor parts. The scalar 
part has the same F dependence [l^l as aj) (F) , and cor- 
rects the latter by approximately 1%, while the second 
rank tensor parts have an F and M dependence, which 
can be parametrized in terms of a third order tensor po- 
larizability aj (6'S'i/2,i^) via 
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FIG. 2: Contribution of off-diagonal hyperfine matrix ele- 
ments to the third order calculation of the Stark effect. 



where the dependence on the angle B between the elec- 
tric field and the quantization axis is given by j(&) — 
3cos^ 6 — 1. The M-dependence in Eq. Q is responsible 
for lifting the Zeeman degeneracy in the ground state hy- 
perfine levels, but gives also a correction to the shift of 
the clock transition frequency which itself is dominated 
by Eq. Q. The third order Stark effect of the two hy- 
perfine levels can then be parametrized in terms of the 
third order polarizabilities by 



a^^\F,M)^4^\F)+a^^\Ff-M^-^^m. 



21(21+1) 



(4) 



In cesium the explicit F-dependence of Eq. Q) yields [T^, 
[13, for 6*^0, 



a(-^^(4,M) = ao + iai + a2) 
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where ai is the contribution of the tensor part of the 
dipole-dipole hyperfine interaction, and 02 the contribu- 
tion of the electric quadrupole interaction. The Fermi- 
contact interaction provides the dominant contribution 
to flo which also has a small contribution from the scalar 
part of the dipole-dipole interaction. Equations Q bear 



a close resemblance to the expressions obtained by San- 
dars f^, except for the negative sign of the ai term in 
Eq. l)5bp . which is positive in Sandars work. We have 
confirmed the sign of our expression in a recent experi- 
ment 11] . That sign error, which seems to have remained 
unnoticed for almost 40 years, has no consequence for the 
tensor polarizability of the F = A state, which is the only 
one measured to date. It affects the static Stark shift 
of the clock transition at a level slightly below today's 
experimental sensitivity. 

Sandars' equations were evaluated by p and who 
considered only diagonal matrix elements of Hhf for the 
QSi/2 and the 6Pj states (Fig.P). They further neglected 
the fine structure energy splitting in the denominators of 

Eq. (Q) and assumed the relation Agp^^, = bAgp^^^ for the FIG. 4: The static Stark shift of the clock transition fre- 

The squares (■) represent experimental values of 
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hyperfine constants, valid for one-electron atoms, while 
for Cs the corresponding ratio of experimental values 
is 5.8. Those assumptions yielded the value (f) in Fig. 
which is in disagreement with the experimental results. 
We have reevaluated |12J | their result by dropping the sim- 
plifying assumptions and by using recent experimental 
values of the reduced matrix elements (6S'i/2 \\d\\ QPj) 
from [lj|. As a result, the discrepancy becomes even 
larger [point (f) in Fig. |3], and does not change signif- 
icantly when the perturbation sum is extended to nPj 
states with n > 6. 

The first term in Eq. QJ is not restricted to diagonal 
matrix elements of the hyperfine interaction. We have 
therefore extended the numerical evaluation of Eq. ^ by 
including off-diagonal terms. Figure [21 gives a schematic 
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FIG. 3: The third order tensor polarizability of the F=4 Cs 
ground state. The filled squares (■) represent experimental 
values of Carrico et al. QCe')! Gould et al. |a|(b), Ospelkaus 
et al. llOlK c'l and Ulzega et al. [01(d). The empty square 
(n) (e) represents a weighted average of (a), (b), (c) and (d). 
The dots (•) represent the theoretical value from [3(f), and 
our re-evaluation of the latter value after dropping simplifying 
assumptions (f). The dotted horizontal line is the result of 
the present work with its uncertainty (shaded band). 



quency. 

Haun et al. [ 3 (a ) , Mowat |l5|l (c'). Simon et al. [l3(e), and 
Godone et al. [13(g)- The dots (•) represent theoretical val- 
ues of Feichtner et al. 1171 (bl. Lee et al. [laKd'l. and Micalizio 
et al. T9l|(f). The circle (o) (b',f') represents the rescaled val- 
ues of (b,f) as explained in the text. The error bar of point 
(e) is smaller than the symbol size. The dotted horizontal line 
shows our result with its uncertainty (shaded band). 



overview of all possible configurations which include off- 
diagonal hyperfine matrix elements and which are com- 
patible with the hyperfine and Stark operator selection 
rules. It is interesting to note that the diagrams 1 and 
2 were already considered by Feichtner et al. [TJ in 
their calculation of the clock transition Stark shift. For 
unknown reasons, the off-diagonal terms were never in- 
cluded in the calculation of the tensor polarizability. 

We have evaluated all the diagrams shown in Fig. [21 
and in particular the diagrams 3, 4, and 5, which, to 
our knowledge, were never considered before. As noted 
in the figure, all diagrams lead to F and M dependent 
energy shifts, except for diagram 1, which gives an M in- 
dependent shift. Thus all diagrams contribute, together 
with the diagonal contributions A and B of Fig. ^ to 
the Stark shift of the clock transition, while only the dia- 
grams B and 2~-5 contribute to the tensor polarizability. 
Moreover, the relative importance of the diagrams for 
the two effects of interest is quite different. In the case 
of the clock shift, we find that 90% (95%) of the total 
contribution (n=6~18) comes from the diagrams A and 
1 evaluated for n = 6,7 [n = 6,7,8). In this case the 
contributing (electric dipole and hyperfine) matrix ele- 
ments are directly or indirectly given by experimentally 
measured quantities. The diagonal hyperfine matrix ele- 
ments are proportional to the measured hyperfine split- 
tings, while the off-diagonal hyperfine matrix elements 
between S states of different principal quantum numbers 
n can be expressed in terms of the geometrical averages 
of the hyperfine splittings of the coupled states, a re- 
lation which holds at a level of 10~^ ^^. This gives 
us a high level of confidence in our value of the clock 
shift rate. The off-diagonal matrix elements of the dia- 
grams 2-5 cannot be traced back to experimental quan- 



titles. We have calculated these matrix elements using 
wave functions obtained from the Schrodinger equation 
for a Thomas-Fermi potential, with corrections includ- 
ing dipolar and quadrupolar core polarization as well as 
spin-orbit interaction with a relativistic correction factor 
following |22|. We included nSj, nPj, and nDj states 
up to n = 18 in the perturbation sum and obtained 



Az^oo/IE^ = -2.06(1) Hz/(kV/cm)2 , 
for the shift of the clock transition frequency, and 



(6) 



.(3) 



a''2'{F ^ I±J)^ =f3. 72(25) x 10"^Hz/(kV/cm)^ , (7) 

for the tensor polarizability. The results are shown as 
dotted hues in Figs. 13 and ^together with previous the- 
oretical and experimental results. The uncertainty of 
our calculated values is indicated by the shaded bands. 
The relative uncertainty of the clock shift is significantly 
smaller than that of the tensor polarizability due to the 
use of experimental values, with relatively small uncer- 
tainties. For the contributions which we calculated ex- 
plicitly from the Schrodinger wave functions we use the 
accuracy (4%-8%) with which these wave functions re- 
produce experimental dipole matrix elements and hyper- 
fine splittings to estimate the precision of our results. 
More details will be given in a forthcoming publica- 
tion [IJ. We can claim that the present calculation of 

(3) 

a2 yields a good agreement with all experimental data. 
The situation with the Stark shift of the clock transi- 
tion frequency is less clear as there are disagreeing exper- 
imental values. While the experimental results Q, [l5l| . 
and [l^ [(a), (c), and (e) in Fig. ^ are supported by 



the theoretical value of [131 (d), our present result is in 
excellent agreement with the recent experimental value 
of Godone et al. [13 (g) and with the calculation of Mi- 
calizio et al. lff|(f). The shift was also calculated by 
Feichtner et al. 17] in an approximation using hydro- 
genic wave functions, neglecting spin-orbit interactions, 
and considering only the diagrams A, B, 1, and 2 [point 
(b) in Fig. ^ . With these approximations the scalar po- 
larizability ao can be factored out of their final result. 
We have rescaled point (b) in Fig^Ejusing more precise 
(consistent) values of ag given in Q,H,I3, yielding point 
(b') which is then consistent with the present result. We 
also rescaled point (f) in the same manner, yielding (f ') 
which cannot be distinguished from (b'). 

We conclude by recalling the relevance of the latter re- 
sult for primary frequency standards. One of the leading 
systematic shifts of the cesium clock frequency is due to 
the interaction of the atoms with the blackbody radia- 
tion (BBR) field. It was shown [221 that the dynamic 
BBR shift can be parametrized in terms of the static 
Stark shift investigated here. The correction factor for 
the BBR shift commonly used is based on the precise ex- 
perimental value of Simon et al. JI3 [point (e) in Fig. 21, 
whose difference with the present result is 21 (53) times 
the corresponding theoretical (experimental) uncertainty. 
In order to shine more light on this important issue we are 
preparing an alternative new experiment for measuring 
the static Stark shift of the clock transition frequency. 
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